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Abstract 

We build on our detailed analysis of time-series observations of the globular cluster M 68 to investigate the irregular 
pulsational behaviour of four of the RR Lyrae stars in this cluster. M 68 is one of only two globular clusters in which 
mode-switching of RR Lyrae stars has previously been reported, and we discuss one additional case, as well as a case 
of irregular behaviour, and we briefly revisit the two previously reported cases with a homogeneous analysis. We find 
that in 2013, V45 was pulsating in the first-overtone mode only, despite being previously reported as a double-mode 
(fundamental and first overtone) pulsator in 1994, and that the amplitude of the fundamental mode in V7 is increasing 
with time. We also suggest that V21 might not have switched pulsation modes as previously reported, although the 


first overtone seems to be becoming less dominant. Finally, 
lost a pulsation mode between 1950 and 1986. 
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1. Introduction 

RRL stars have been the subject of intense study for more 
than a century, meaning that for many of these objects, 
time-series photometric data are available with long base¬ 
lines. This allows for the study of any changes in the pul¬ 
sational behaviour, including period changes, Blazhko-type 
modulation, and mode switching. This is particularly true 
for RRL stars in globular clusters, as clusters are an ideal 
environment for stellar population studies. Indeed, observ¬ 
ing campaigns to monitor globular clusters are still ongoing 
both from the ground (e.g. Kains et al. 2015, 2013; Arellano 
Ferro et al. 2015; Carretta et al. 2014; Kunder et al. 2013; 
Amigo et al. 2013; Figuera Jaimes et al. 2013) and from 
space (e.g. Neeley et al. 2015; Contreras Ramos et al. 2013; 
Montiel & Mighell 2010). 

Only a few cases of RRL stars switching their pulsa¬ 
tion modes have been reported in the literature. The first 
cases were reported in the cluster that is the subject of 
this paper, M 68, by Clement et al. (1993, hereafter C93). 
Using observations by van Agt & Oosterhoff (1959, here¬ 
after V059) from 1950 and Rosino & Pietra (1954, here¬ 
after R54), made between 1951 and 1953, Clement (1990) 
found evidence that V33 was pulsating as a RROl star^. 

^ In this paper we use the notation of Alcock et al. (2000) to 
refer to single-mode fundamental and first-overtone pulsators as 


our analysis of available archival data confirms that V33 


However, in observations from 1986-1991 taken and anal¬ 
ysed by C93, V33 was found to have changed into a RRl 
variable. In later years, Brocato et al. (1994, hereafter B94), 
Walker (1994, hereafter W94) and Kains et al. (2015, here¬ 
after K15), found that V33 continued to pulsate in the first 
overtone (FO) mode during observations from 1989, 1993, 
and 2013, respectively. C93 also found that V21 was pul¬ 
sating in both the fundamental (F) and FO modes with 
approximately equal amplitudes over three years from 1986 
to 1988, and that from 1989 to 1991 it was pulsating in the 
FO mode only, with no evidence for F mode pulsations. 

The variable V79 in M 3 was the next case of mode¬ 
switching to be found. Clement & Shelton (1997) discovered 
that it had changed from RRO in 1920-1926 (Larink 1922; 
Muller 1933; Greenstein 1935; Nemec & Clement 1989) to 
RROl in 1996, but it was then found in the 2008 observa¬ 
tions of Goranskij et al. (2010) to have reverted back to 
single-mode pulsation in the F mode (RRO). Also in M 3, 
Corwin et al. (1999) found that V166 had switched from 
F to FO pulsation over the course of one year based on 
observing runs in 1992 and 1993, while Clementini et al. 
(2004) found that two double-mode variables in the same 
cluster, V200 and V251, had changed their dominant pulsa- 


RRO and RRl, respectively, and double-mode (fundamental and 
first-overtone) pulsators as RROl. 
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tion mode over the course of the same time interval, using 
the data of Corwin et al. (1999), but remained RROl vari¬ 
ables. 

Recently, Soszyhski et al. (2014a) reported on mode¬ 
switching in OGLE-BLG-RRLYR-12245, an RRL star that 
was first discovered from OGLE-111 observations taking 
place between 2001 and 2006, and subsequently re-observed 
with the start of OGLE-IV observations. Whereas OGLE- 
BLG-RRLYR-12245 was a double-mode RROl pulsator 
in the 2001-2006 observations, the 2010-2013 light curve 
shows clearly that the star was pulsating only in the E 
(RRO) mode. The light curve from Soszyhski et al. (2014a) 
also reveals a slow increase in amplitude in the 2001- 
2006 light curve. Erom the same data, Soszyhski et al. 
(2014b) found two more mode-switching stars: OGLE- 
BLG-RRLYR-07226, which also went from RROl to RRO, 
and OGLE-BLG-RRLYR-13342, which changed from RRO 
to RROl. In the LMG, Poleski (2014) found that OGLE- 
LMC-RRLYR-13308 had switched pulsation modes from 
RROl to RRO, using OGLE and EROS-11 photometry. 
Einally, Drake et al. (2014) found 6 more mode-switching 
RROl stars in the Catalina Surveys of Periodic Variable 
Stars, in this case changing their dominant pulsation mode 
over just a few months, although at least one of these is in 
fact not a case of mode-switching (Soszyhski et al. 2014b). 

Although there is general agreement, backed up by ob¬ 
servational evidence, that RROl stars are associated with 
the region of the horizontal branch (HB) in which stars can 
pulsate either in the P or EO mode (Soszyhski et al. 2009), 
there is no consensus on the exact causes for double-mode 
pulsation (e.g. Smolec 2014; Dziembowski 1993). Studies 
including turbulent convection models were once thought 
to yield stable double-mode pulsation (e.g. Szabo et al. 
2004), but Smolec & Moskalik (2008) showed that this was 
due to incorrect treatment of buoyant forces, and were un¬ 
able to reproduce these results with the correct analysis. 
Furthermore, not all stars located in the either-or region 
are double-mode pulsators, as non-linear effect are respon¬ 
sible for the double-mode pulsation, meaning that the locus 
of RROl stars should be narrower than the either-or region. 
Interestingly, however, this is not supported by the observa¬ 
tions of K15 (see their Fig. 11, or Fig. 1 in this paper), who 
found that the RROl stars in M 68 span most of the colour 
range of the instability strip. Indeed, Arellano Ferro et al. 
(2015) argued that the separation of the RRO and RRl pop¬ 
ulations (shown in Fig. 1 as a vertical dotted line), corre¬ 
sponds to the first overtone region’s red edge (FORE), and 
therefore any double-mode star that is redder than this em¬ 
pirical border cannot be in the either-or region. Increasing 
the sample of detected mode-switching RRL stars is an im¬ 
portant step towards understanding whether this behaviour 
is typical for stars evolving along the HB. This would also 
allow us to constrain the changes in internal structure that 
these stars are going through when moving across the in¬ 
stability strip, which could in turn lead to clearer explana¬ 
tions for other pulsational irregularities such as the Blazhko 
effect. Here, we report at least one additional case of mode¬ 
switching in M 68, placing this cluster on a par with M 3 
with regards to the number of known cases. 

M 68 (NGC 4590, C1236-264 in the lAU nomenclature; 
a = 12'‘39'"27.98*,5 = -26°44'38.6" at J2000.0) is one of 
the most metal-poor globular clusters in the Galaxy, and 
contains 50 known variable stars. Recently, K15 studied the 
variables in this cluster in detail, and used Fourier decom¬ 


position of the RR Lyrae (RRL) and period-luminosity rela¬ 
tions for the SX Phoenicis stars to estimate the metallicity 
and distance of the cluster. The 42 RRL variables in the 
cluster include 12 double-mode RRL (RROl) stars, which 
are particularly interesting, because the double-mode pul¬ 
sation allows us to measure the mass of these objects (e.g. 
Petersen 1973; Popielski et al. 2000; Marconi et al. 2015), 
as well as to compare the mass-metallicity distribution of 
field and cluster RRL stars (e.g. Bragaglia et al. 2001). 

In this research note, we report irregularities in the be¬ 
haviour of some of the variables using archival data sets 
dating back over 65 years. We provide a brief summary 
of our observations in Sec. 2, and discuss one new case of 
mode switching in M 68 in Sec. 3, and irregular behaviour 
of one of the double-mode pulsators in Sec. 4. In Sec. 5 we 
revisit the two previously reported cases of mode switching 
in M 68, and we briefly summarise our findings in Sec. 6. 
The light curve data used in this paper are available for 
download via the GDS^, and are linked via the relevant 
publications on the NASA ADS system. 

2. Observations and reductions 

Images were obtained in V (236) and I (219) using 
GGD cameras on the Im LGOGT telescopes at the South 
African Astronomical Observatory (SAAO) in Sutherland, 
South Africa, as well as the Gerro Tololo Inter-American 
Observatory (CTIO) in Ghile. The observations covered 74 
days from 8th March 2013 until 20th May 2013. The im¬ 
ages were reduced using the DanDIA difference image analy¬ 
sis software (DIA; Bramich et al. 2013; Bramich 2008), and 
light curves were obtained for all stars detected in reference 
images built from the best-seeing images in each filter. For 
further details on the observations, reductions, and calibra¬ 
tion of the photometry, the interested reader is referred to 
K15. 

3. \/45: A new case of mode-switching in M 68 

V45 was hrst detected by W94, who classified it as a 
double-mode RROl star, with periods Pq = 0.5239 d and 
Pi = 0.39087 d, giving a period ratio of 0.746. However, 
K15 found no evidence of a second pulsation mode in their 
recent data (Table 1). Comparing the light curves of W94 
with those of K15 (Fig. 2) suggests indeed that the double¬ 
mode pulsation that was present in the 1993 observations 
essentially disappeared between 1993 and 2013, with the 
most recent light curve looking like that of an RRl star 
with period of 0.3908187 d. The position of V45 on the HB 
was significantly different in the data of K15 compared to 
the observations of W94. Its location near several bright 
stars (see Fig. 10 of K15) means that it is likely that W94 
over-estimated the brightness of V45 because of blending, 
which had little effect on the photometry of K15 thanks 
to their use of difference image analysis (e.g. Kains et al. 
2012; Bramich et al. 2012). We note here that K15 under¬ 
estimated the average brightness of V45 in V by 0.04 mag, 
and therefore its H — / colour. Here we revise this value to 
find ly) = 15.68 ±0.01 mag, and V — I = 0.51 ±0.01 mag, 
placing it among the RRl population. 

Performing a Fourier analysis of the data of W94, we 
found an amplitude for the F mode pulsation of 0.045 mag 

^ http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/578/A128 
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Data 

Observed 

(year) 

Ab 

(mag) 

Av 

(mag) 

Aj 

(mag) 

Po 

(d) 

Pi 

(d) 

ks 

ky 

ki 

Type 

W94 

1993 

0.74 

-0.47 

-0.45 

0.5076(2) 

0.38450(4) 

4.7(9) 

3.0(6) 


RROl 

K15 

2013 

- 

0.52 

0.30 

- 

0.39082(1) 

- 

- 

- 

RRl 


Table 1: Same as Table 2, but for V45. 
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Figure 1: The HB of M 68, with RRL stars overplotted. 
The entire RRL population is shown in grey, with RRO as 
filled circles, RRl as filled triangles, and RROl as inverted 
open triangles. The four stars discussed in this paper are 
shown as red filled squares (K15 observations) and open 
green circles (W94), with dashed lines linking their K15 
and W94 positions. The vertical dotted line indicates the 
separation between the loci of RRO and RRl stars. 
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Figured: Phased light curves (with the FO periods given in 
Table 2) for V7 showing the data of V059 (unfiltered pho¬ 
tographic plates, top), R54 (unfiltered photographic plates, 
middle left), C93 (B-band, middle right), W94 (fo, bottom 
left), and K15 {V, bottom right). The data of B94 are not 
plotted because they have poor phase coverage. 


Figure 2: fo—band light curves for V45 showing the data of 
W94 (left) and K15 (right), phased with the FO periods 
listed in Table 1. 


in V and 0.022 mag in I. The residual Fourier noise of the 
K15 light curve, after prewhitening the FO pulsation, is 
0.008 mag in V, and 0.006 mag in I (see Table 3), mean¬ 
ing that if present in 2013, the F-mode pulsation would 
have been detected in both filters, unless the F mode had 
dropped below detection limits. 

4. Amplitude ratio variations in V7 

V7 was discovered by Shapley (1919), and was first identi¬ 
fied as a double-mode pulsator by C93. K15 calculated F 
and FO periods of Pq = 0.520186 d and Pi = 0.3879608 
d, with a period ratio of 0.746, typical of RROl pulsators 
(e.g Poleski 2014), and similar to values found by C93 and 
W94. Its position on the HB, along with the other three 
stars discussed in this paper, is shown in Fig. 1. 

Comparing our data with archival data, we noticed that 
the photographic plate light curve amplitudes of V059 and 


R54 are 0.33 and 0.36 mag, respectively, whereas the am¬ 
plitude in B found by C93 was 0.87 mag, and in V, K15 
found an amplitude of 0.62 mag, similar to the V ampli¬ 
tude of 0.63 mag found by W94. This is surprising, because 
amplitudes measured from photographic plates should be 
larger than in the redder V hlter. This suggests that the am¬ 
plitude of one of the pulsation modes has grown in the time 
between the observations of V059 and R54, and those of 
C93, W94, and K15. The photographic light curves taken in 
the 1950s do not show as large scatter due to double-mode 
pulsation as the more recent data sets, and R54 classified 
V7 as a RRl star. 

The results of a homogeneous re-analysis of all avail¬ 
able light curves for V7 with the period analysis software 
Period04 (Lenz & Breger 2005) are given in Table 2. The 
ratio of the FO to F mode amplitudes, k, is given for each 
data set and each filter. This ratio decreases with time, 
supporting our hypothesis that the F mode is growing over 
time. By comparing the amplitudes of the two modes, we 
note, however, that the amplitude of the F mode for the 
case of K15 is ^^^0.4 mag in V, which is too small for the 
F mode period when compared to the locus of RRO stars 
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Data 

Observed 

(year) 

Ap 

(mag) 

Ab 

(mag) 

Av 

(mag) 

Al 

(mag) 

Po 

(d) 

Pi 

(d) 

kp 

ks 

ky 

ki 

Type 

V059 

1950 

0.33 

— 

— 

- 

0.5174(1) 

0.3881(6) 

3.0(5) 

— 

— 

— 

RROl 

R54 

1951-1953 

0.36 

— 

— 

- 

0.54118(4) 

0.387951(5) 

5.1(9) 

— 

— 

— 

RROl? 

C93 

1986-1991 

- 

0.87 

- 

- 

0.519742(8) 

0.387958(2) 

- 

3.0(3) 

- 

- 

RROl 

W94 

1993 

- 

0.87 

0.63 

0.53 

0.51929(8) 

0.38801(2) 

- 

3.0(3) 

2.5(2) 

1.4(1) 

RROl 

K15 

2013 

- 

- 

0.62 

0.41 

0.52028(8) 

0.38812(2) 

- 

- 

1.8(1) 

1.7(1) 

RROl 


Table 2: Light curve amplitudes (peak-to-peak for the filtered light curves, Fourier amplitudes for the photographic data) 
of the phased light curves and periods of the double-mode variable V7 in the different available data sets. All periods 
have been recalculated with the software Period04, and amplitudes were also re-measured from the light curves. The 
ratio of the FO to F amplitudes, k, is also given for each filter, and p denotes photographic data. We do not include 
B94 data because of poor phase coverage. The number is parentheses next to each value is the uncertainty on the last 
decimal point. 


on the period-amplitude diagram of M 68 (Fig. 17 of K15). 
This indicates that the amplitudes of each mode in double¬ 
mode pulsators are smaller than their single-mode coun¬ 
terparts, which is expected if non-resonant mode coupling 
takes place (e.g. Smolec 2014), as is expected to be the 
case in a majority of double-mode stars. In any case, even 
if this indicates that the F mode is weak, it is still becoming 
stronger relative to the FO mode. 

In /, the ratio increased from W94 to K15, but we be¬ 
lieve kj to be underestimated in the data of W94 due to 
a poor light curve. Although we do detect evidence for F 
mode pulsation in the R54 data, we find a period for this 
mode of 0.541177 d, which is much longer than the funda¬ 
mental period derived from any of the other data sets. This 
could indicate that a mode switch was about to happen, or 
had recently happened, in a similar manner to V79 in M 
3 (Goranskij et al. 2010), but we lack the observations to 
verify this possibility. Furthermore, the F mode detection 
is weak, with a ratio k of 5.16, leaving open the possibility 
that V7 might have been pulsating as a RRl star in the 
observations of R54. 

We note that V7 has the shortest periods (both Pq and 
Pi) of all RROl stars in M 68 (K15). The periods derived 
for the pulsation inodes of V7 also seem to vary significantly 
with time, particularly the F mode period. 
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5. Previously reported cases 

5.1. V21 

As mentioned in Sec. 1, V21 was found by C93 to have 
switched from being an RROl pulsator between 1986 and 
1988 to an RRl star between 1989 and 1991, meaning that 
the transition must have been abrupt. However, both W94 
and K15 found that V21 was once again a clear double¬ 
mode RROl pulsator, meaning that another transition must 
then have occurred between 1991 and 1993. We calculated 
Fourier fits using Period04 in order to determine the am¬ 
plitude of each mode in the data sets of C93, W94, and K15. 
We also considered the 1986-1988 and 1989-1991 parts of 
the C93 light curve separately, to which we refer as C93A 
and C93B respectively. 

Although our analysis does confirm that two modes are 
detected in C93A, and only one in C93B (Table 4), we 
suggest that the poor quality of the C93 light curve (Fig. 4) 
does not actually allow a strong conclusion to be made 
as to the pulsational behaviour of V21. The evidence for 
an abrupt mode switch from RROl to RRl and then back 


Figured: Light curves for V21 showing the data of V059 
(unfiltered photographic, top left), C93 (P, top right, C93A 
and C93B are shown separately in the middle panels), W94 
{V, bottom left), and K15 (R, bottom right). The light 
curves are phased with the FO periods listed in Table 4. 


to RROl over the course of 5 years is therefore weak. In 
addition to this, analysis of the C93A data suggests that 
the F mode was dominant over 1986-1988, which is unusual 
for an RROl pulsator (Smolec et al. 2015), and the ratio of 
periods of the detected modes is 0.775, significantly higher 
than the canonical value of ^0.746. 

However, there is evidence of interesting behaviour 
in the V21 light curves. We re-analysed the light curves 
of V059 and were unable to find signs of the F mode 
pulsation, although the light curve is likely too sparse to 
draw conclusions as to the mode content of V21 in 1950. 
We also found that the amplitude ratio k decreased in both 
V and I between the W94 and K15 data. Since very similar 
light curve amplitudes were found in those two studies in 
both V and /, along with very similar periods (Table 4), 
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V7 

V21 

V33 

V45 

Data 

Filter 

FO 

F 

FO 

F 

FO F 

FO 

F 

V059 

P 

0.011 

0.012 

0.038 

- 

— — 

- 

— 

R54 

P 

0.013 

0.012 

— 

— 

— — 

— 

— 

C93A 

B 

— 

— 

0.037 

0.033 

— — 

— 

— 

C93B 

B 

— 

— 

0.049 

— 

— — 

— 

— 

C93 

B 

0.012 

0.011 

0.040 

0.042 

0.012 - 

— 

— 

W94 

B 

0.014 

0.012 

0.018 

0.018 

0.009 - 

0.025 

0.024 

W94 

V 

0.010 

0.009 

0.014 

0.014 

0.007 - 

0.017 

0.014 

W94 

1 

0.013 

0.012 

0.008 

0.009 

0.005 - 

0.016 

0.016 

K15 

V 

0.015 

0.015 

0.016 

0.016 

0.009 - 

0.008 

— 

K15 

I 

0.010 

0.010 

0.010 

0.010 

0.004 - 

0.006 

- 


Table 3: Fourier noise after prewhitening, in magnitudes, at 
this paper, for each filter and each data set. 



'll 

Figure 5: Light curves for V33 showing the data of V059 
(unfiltered photographic, top left), C93 (B, top right), W94 
(F, bottom left), and K15 (F, bottom right). The light 
curves are phased with the FO periods listed in Table 5. 


this change is likely to be real, suggesting that the relative 
amplitude of the FO to the F mode has decreased between 
1993 and 2013. 


5.2. V33 

We re-analysed the data of V059 for this star using 
Period04, and confirm that there are signs of double-mode 
pulsation in that data set, as first suggested by Clement 
(1990). We also find a FO period for the V059 data that 
is significantly shorter than those subsequently found by 
C93, W94 and K15, when the star was a RRl pulsator 
(Table 5). Therefore, we confirm that V33 switched from 
being a RROl to RRl star between 1950 and 1986. We also 
note that the FO period is becoming longer with time. 

6. Conclusions 

Further comparative analysis of the data of K15 with 
archival data revealed interesting behaviour in several RRL 
stars. The F mode in V45 has either disappeared or become 
very weak, meaning that it now pulsates like a RRl star, 
with the switch having occurred between 1993 and 2013. 


the F and FO frequencies for the four variables discussed in 


We find that the F mode in V7 is growing, and that a 
mode-switching event might have taken place in the 1950s, 
associated with a large change in the F pulsation period, 
but the evidence is too weak to make strong conclusions 
on this for V7. The F mode of V21 has grown with time, 
although again the claim of mode switching is weak, as it 
relies solely on a light curve that suffers from large scatter. 
We also find evidence to support the claim of C93 that V33 
switched modes from RROl to RRl between 1950 and 1986. 
Interestingly, the double-mode stars in M 68 are distributed 
over a large part of the instability strip, which does not sup¬ 
port the current understanding of double-mode pulsation. 
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Table 4: Same as Table 2, but for V21. 
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Table 5: Same as Table 2, but for V33. No amplitude is given for the V059 data as they are not on a magnitude scale. 
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